Optimized coupling of cold atoms into a fiber using a blue-detuned hollow-beam funnel 
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We theoretically investigate the process of coupling cold atoms into the core of a hollow-core 
photonic-crystal optical fiber using a blue-detuned Laguerre-Gaussian beam. In contrast to the use 
of a red-detuned Gaussian beam to couple the atoms, the blue-detuned hollow-beam can confine 
cold atoms to the darkest regions of the beam thereby minimizing shifts in the internal states and 
making the guide highly robust to heating effects. This single optical beam is used as both a funnel 
and guide to maximize the number of atoms into the fiber. In the proposed experiment, Rb atoms 
are loaded into a magneto-optical trap (MOT) above a vertically-oriented optical fiber. We observe 
a gravito-optical trapping effect for atoms with high orbital momentum around the trap axis, which 
prevents atoms from coupling to the fiber: these atoms lack the kinetic energy to escape the potential 
and are thus trapped in the laser funnel indefinitely. We find that by reducing the dipolar force 
to the point at which the trapping effect just vanishes, it is possible to optimize the coupling of 
atoms into the fiber. Our simulations predict that by using a low-power (2.5 mW) and far-detuned 
(300 GHz) Laguerre-Gaussian beam with a 20-/im radius core hollow-fiber it is possible to couple 
11% of the atoms from a MOT 9 mm away from the fiber. When MOT is positioned further away, 
coupling efficiencies over 50% can be achieved with larger core fibers. 

PACS numbers: 37.10.Gh 



I. INTRODUCTION 

The field of atom optics has produced matter-wave 
equivalents for almost every conventional optical compo- 
nent [1-3]. However, two key elements are still missing: 
(i) a flexible atomic waveguide that acts as the equivalent 
of an optical fiber, together with (ii) the tools necessary 
for coupling free-space atoms into these flexible waveg- 
uides. Efficient guiding in optical fibers would allow for 
delivery of atoms over long distances and arbitrary paths 
that might enable new applications in, e.g., atom inter- 
fcrometry [4, 5], gravity sensors [6], ultra-precise atom 
implantation [7, 8], and atomic lithography [9]. 

In this paper, we develop a solution for both elements 
with a single optical beam. We consider the particular 
situation in which we use a repulsive potential to guide 
the cold atoms. This approach is interesting because the 
atoms can be guided in a minimally perturbating loca- 
tion with minimal unwanted heating. Unfortunately, as 
we will show, this approach has potentially much lower 
coupling efiiciency than the opposite approach using an 
attractive potential. However, we will identify the con- 
ditions in which this repulsive guiding can provide the 
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same order of efficiency as the attractive potential. We 
also show a rather unexpected conclusion that there is 
an optimal potential strength, under a given set of initial 
conditions, to maximize the atomic coupling. In other 
words, the deepest possible potential does not provide the 
best coupling of free-space atoms into the optical fiber. 

Following the first theoretical proposals [10, 11], room- 
temperature atom guidance in glass capillaries was 
quickly demonstrated using red-detuned Gaussian beams 
[12] as well as blue-detuned evanescent waves [13, 14]. In 
contrast to guiding with red-detuned beams, which draw 
atoms into the high intensity areas of the beam, blue- 
detuned guides present the inherent advantage of min- 
imizing the heating of atoms due to spontaneous emis- 
sion and reduce perturbation of the energy levels by the 
light potential (Stark shifting) [11]. Shortly after these 
first demonstrations, evanescent guiding of laser-cooled 
atoms using beams of He [15] and Rb [16] was achieved. 
Further progress was limited because optical guidance in 
capillaries is both lossy [17] and highly multimode [15] 
which leads to high atomic losses and short guidance 
lengths (^cm). One of these models [18] also proposed 
the use the diffracting field from the cladding to cou- 
ple the atoms, although the low power in the evanescent 
field required a high-power near-detuned beam with its 
associated spontaneous emission losses. 

Development of hollow-core photonic-crystal fibers 
(HC-PCF) [19, 20] paved the way for efficient guidance of 



atoms because of their low loss and single mode optical 
guiding. The first demonstration of room-temperature 
atoms guiding through a HC-PCF reported a guid- 
ing efficiency of up to 70% through a 6-cm-long fiber 
[21]. Following soon after, two experiments demon- 
strated the loading of cold atoms into HC-PCF using 
red-detuned Gaussian beams [22. 23]. Very recently, cold 
Rb atoms were guided through a 8.8-cm-long HC-PCF by 
a 2.3 Watt far-detuned Gaussian beam and a peak flux 
of 10^ atom/s was observed [24]. 

Two approaches were recently developed to improve 
the coupling efficiency of cold atoms into small-core HC- 
PCF. Bajcsy et al. used a quadrupole magnetic funnel 
to guide the MOT from its launch distance to within 
reach of the Gaussian beam funnel [23]. Soon after, they 
improved the atom coupling efficiency by a factor of 6 
using a collimated hollow-beam and a repulsive sheet 
beam. The system effectively acts as an "optical eleva- 
tor" that lowers the MOT within 1 mm of the fiber input 
[25]. When the sheet beam is turned off, gravity and a 
red-detuned Gaussian beam funnels the atoms into the 
core. 

For blue-detuned atomic guidance, HC-PCF enables 
the use of a hollow Laguerre-Gaussian beam itself to 
guide the atoms instead of just the evanescent-field part 
that leaks into the core. This allows for a deeper trapping 
potential at much lower optical power, while keeping the 
atoms away from influence of the fiber walls [11, 26]. Col- 
limated Laguerre-Gaussian (LG) beams have also been 
used to guide cold atoms from a MOT and it was shown 
that high ^-order beams proved to be the most efficient 
guides as they minimize the potential energy given to 
atoms at the loading stage [27]. However, to our knowl- 
edge, no HC-PCF have been demonstrated to guide high 
£-order LG beams. On the other hand, low-loss guid- 
ance of a LGoi-like mode in a 19-m-long HC-PCF was 
recently demonstrated in [28]. Models of atom guiding 
in collimated blue-detuned hollow-beams in vacuum have 
also been developed and validated with experimental re- 
sults [27, 29; 30], although to our knowledge no one has 
yet considered the cold atoms behavior in the diffracted 
field from the core of the fiber. 

In this paper, we present the solution to the best pa- 
rameters for optimal atomic coupling. In Sec. II, the 
physical equations and 3D Monte Carlo numerical model 
of our proposed experiment are presented. Sec. Ill deals 
with the dynamics of the atoms in the funnel and the 
gravito-optical bottle trapping effect. The system's atom 
coupling efficiency against experimental parameters is 
studied in Sec. IV. 
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FIG. 1. (Color online) Experimental setup for coupling atoms 
from a MOT into a hollow-core photonic-crystal fiber. The 
fiber is mounted vertically in the vacuum chamber. A single 
blue-detuned collimated Laguerre-Gaussian beam is coupled 
into the fiber core at the bottom of the fiber. At the top end, 
the beam difitracts into free space acting as a funnel. When 
the MOT is released, from distance d of the fiber output, the 
combined action of gravity and the hollow-beam funnels the 
atoms into the core and guides them down the fiber. The 
schematic is not drawn to scale. 



its axis is parallel to gravity. A blue-detuned hollow- 
beam diffracting from the core acts as a funnel to both 
couple and guide atoms released from a MOT. 

In this paper, the term "coupling" is used to describe 
the process of optically funnelling the atoms into the fiber 
core. The term "guiding" is exclusively used to describe 
the optical control of the motion of atoms when they are 
inside the fiber core. 



A. Forces in the optical funnel 



II. 



ATOM COUPLING MODEL 



We propose, in Fig. 1, an experimental arrangement in 
which a far-detuned low-order Laguerre-Gaussian beam 
is coupled and guided in a hollow-core photonic-crystal 
optical fiber [19]. The fiber is aligned vertically, so that 



The rotational symmetry allows the problem to be 
analysed using cylindrical coordinates: r being the dis- 
tance to the fiber core central axis and z the axial dis- 
tance from the fiber. The extremity of the fiber fac- 
ing the MOT defines the axial origin (z = 0). We will 
use this form to describe the intensity of a single-ringed 



Lagucrrc-Gaussian mode (radial order p = 0, azimuthal 
order ^^0) [31]: 
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where P is the 



optical power and w{z) = 
woa/1 + (z/zfl)^ is the radius of the lowest order 
beam propagating in free space. The beam waist radius, 
inside the fiber, is wq and zn = ttwq/A is the Rayleigh 
range. When the laser frequency is far-detuned (>1000 
r, >40GHz for Rb atoms), only the "off-resonant" 
dipolar force is effective and the trapping potential is 
expressed in the two- level approximation [32]: 
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where the laser detuning S = ujiaser — wq, F is the nat- 
ural linewidth of the cooling transition, and Is is the 
saturation intensity parameter. In this paper, we will 
specify the energies of atoms in terms of an equivalent 
thermal energy for case of comparison between energy 
scales. The dipolar force is equal to the negative gradi- 
ent of the potential C/(r, z). In the interest of readability, 
the vector force will be separated into its axial {F^} and 
radial {Fr) components. Explicit mention of the vari- 
able's dependence in r,z will also be omitted. Deriving 
the trapping potential [Eq. (2)] in the context of the LG 
intensity [Eq. (1)], in r and 2, we obtain: 
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As an example, when comparing the axial and transverse 
components of the dipolar force at the same spatial co- 
ordinates in an LGqi beam with a l/e^ waist diameter of 
40 /im (diffraction half-angle of 1.6°), we find that Fz is 
at least 3 orders of magnitude weaker than Fj. . Nonethe- 
less, Fz can be many times stronger than gravity, even for 
a realistic trap potential depth of 0.5 mK. This results in 
Fz having a surprisingly strong influence on the dynam- 
ics of the atoms in the trap. It can be seen from Eq. (2) 
that for large dctunings, the potential is merely linearly 
dependent upon the ratio of the optical power over the 
detuning. From now on, we will refer to this parameter 
as the light force parameter k = P/5. 

The radial force confines atoms within the funnel as 
gravity pushes them downwards toward the fiber core. 
Fig. 2 shows the radial acceleration on *^Rb resulting 
from a dipolar radial force under typical conditions at 
three axial distances from the end of the fiber. The force 
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FIG. 2. (Color online) Radial component of the blue-detuned 
dipolar force [see Eq. (3)] at three distances from the fiber end: 
mm (blue dashed), 0.4 mm (red solid) and 2.4 mm (black 
dash-dotted). The force is represented as ®^Rb acceleration 
in units of standard gravity. Positive values means the force 
vector is oriented towards increasing values of r (right), while 
negative values indicate the force vector points to decreasing 
values of r (left). This calculation used a 4mW LGoi beam, 
blue-detuned to 300 GHz from the ^^Rh D2 cooling transi- 
tion {5^Si/2iF = 3) ^ 5^P3/2(i^ = 4)) around 780.24 nm 
and diffracting from a 40-/im-diameter fiber core. The tran- 
sition natural linewidth is r/27r = 5.89 MHz and saturation 
intensity is Isat = 1.63niW/cm^ [33]. 



has the form of a trapping potential, which has its max- 
imum at the output of the fiber and decreases with dis- 
tance proportional to the intensity. An important con- 
sequence of the intensity distribution of a blue-detuned 
hollow-beam is its limited confinement radius. Outside 
the maximum intensity radius, the decreasing gradient 
produces a force that repels atoms away from the axis of 
the beam. 



The axial force, i^^, is a manifestation of the gradient 
of intensity produced by a diverging beam. Therefore, it 
is null for collimated beams, e.g. inside the fiber. The 
axial force is always directed away from the end of the 
fiber for a blue-detuned beam regardless of the optical 
beam's propagation direction. This force is seen to rise 
sharply when the atom conies within 1.5 zn of the end 
of the fiber. Under typical coupling conditions, the axial 
force will be many times the gravitational force when 
close to the radius of maximum intensity. As an example. 
Fig. 3 shows the axial acceleration on a ®^Rb atom for 
three axial distances, z: 0.02, 0.4 and 2.4mm. The axial 
acceleration reaches a maximum at 0.4 mm. At distances 
of 0.02 and 2.4 mm, the axial force maximum is equal to 
the gravity. 
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FIG. 3. (Color online) Axial component of the blue-detuned 
dipolar force [see Eq. (4)] at three distances from the fiber end: 
0.02mm (blue dotted), 0.4 mm (red solid) and 2.4 mm (black 
dashed). The force is represented as ^^Rb acceleration in units 
of standard gravity. Positive values means repulsive while 
negative is an attractive force. The same physical parameters 
have been used as in Fig. 2. 



B. Monte Carlo simulation model 
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FIG. 4. (Color online) (a) A characteristic orbital trajectory 
of an atom captured by the confining potential and coupled 
into the fiber, (b) Top- view of the elliptical orbital trajectory 
as the atom is guided toward and inside the fiber. By conser- 
vation of angular momentum, the atom accelerates transver- 
sally as its elliptical movement is compressed into an ever 
smaller orbital. 



A key objective of this research was to identify the 
best coupling efficiency using blue-detuned light with re- 
alistic representations of the initial state. The coupling 
efficiency is defined as the ratio of the number of atoms 
that enter the fiber core over the total number of atoms 
loaded in the MOT. The Monte Carlo method was used 
to determine an initial state for each atom in the MOT. 
The atoms are treated as independent point particles in 
the MOT and we used a Maxwell-Boltzmann distribution 
for their velocities [34] . We selected a temperature of 25 
/iK, which can be achieved with careful tuning of a stan- 
dard MOT [35]. We modelled the density variation as 
Gaussian {p{r) = poexp{—r^/a^)) with the 1/e radius of 
the trap of 70 /xm, which can be experimentally achieved 
with careful tuning of the optical and magnetic forces to 
optimize the trap spring constant [35, 36]. 

We numerically simulate the trajectory of each atom in 
the MOT from the time-of-release. The Velocity- Verlet 
algorithm was selected because of the algorithm's prop- 
erty to conserve the energy of particles modelled in a 
conservative system [37]. At every timestep At, the po- 
sition and velocity of each atom is calculated as follow: 

r{t + At)=r{t)+v{t)At + ^a{t)At'^ (5) 

v{t + At) = v{t) + i [a{t) + a{t + At)] At (6) 

The accelerations a{t) and a{t + At) are both simply the 
sum of the accelerations by gravity and the dipolar force 
Eqs. (3), (4) at their corresponding positions in time r(i) 
and r{t + At). Errors in the simulation can be readily 
estimated by observing changes in the particle energy 
during the simulation. 



III. 



ATOM COUPLING DYNAMICS 



In this section, we undertake a detailed examination of 
the atom trajectories to understand the behaviour of the 
atoms. We distinguish the outcome into three classes: 
(i) atoms coupled into the fiber ("coupled"), (ii) those 
captured in the funnel but which do not enter the fiber 
("trapped"), (iii) and those which immediately escape 
the funnel when released from the MOT ("escaped"). 
In Sec. Ill A, we will analyze atomic motion and then, 
adding the particular details of a MOT released into an 
optical blue-detuned funnel in Sec. IIIB. 



A. Motion of atoms in an optical dipolar potential 

By considering atoms inside the fiber, whether guided 
in a blue-detuned collimated Laguerre-Gaussian beam or 
red-detuned Gaussian beam, their motion will be that of 
a particle trapped in a central, attractive and circular 
potential similar to the case of a test mass moving in 
the gravitational field of a point mass. More generally, 
for a central power-law potential V = r'^ , only potentials 
in a; ~ —1 (gravity) or a; = 2 (harmonic) can sustain 
stable closed elliptical orbitals [38] . The optical confining 
potential follows a sigmoidal-shaped function of r e~^ . 
Therefore, we expect that atoms guided/loaded in the 
optical dipolar potential will not describe stable orbitals, 
which we do find, as shown in Fig. 4. This figure displays 
a characteristic orbital for an atom falling into the funnel. 

If we consider an atom falling into a diffracting poten- 



tial, the atom will spiral into the light funnel, analogous 
to a ball moving on a hyperbolic funnel surface. As the 
diameter of the potential reduces towards the waist of 
the beam, the orbital velocity of the moving particle is 
seen to increase as is expected from conservation of the 
angular momentum. Unlike the ideal circular motion in 
the classic hyperbolic funnel, we see the orbitals describe 
open precessing ellipses as seen in Fig. 4. Similar trajec- 
tories are also observed in collimated beam simulations. 
Depending on the initial conditions, simulated orbitals 
range from almost circular to extreme ellipses (close to a 
purely radial motion). 

In the optical funnel, the axial motion is linked to the 
transverse velocity in a complex fashion. In the absence 
of axial forces in the funnel, i.e. where the diffraction is 
absent or inside the fiber, all atoms accelerate downwards 
with gravity independently of their orbital motion. How- 
ever, in locations where diffraction is significant, the tan- 
gential velocity assumes much greater importance. The 
orbital radius of a high tangential velocity atom needs 
to be larger in order to be balanced by the stronger ra- 
dial forces available at that location. In this case the 
atom will penetrate further into the high intensity areas 
of the beam and thus, the effective vertical acceleration 
delivered by the hollow-beam will be stronger. By this 
means, the optical funnel can couple tangential and axial 
motions of the atom and yield the complex orbits that 
are seen on Fig. 4 as well as the axial velocities illustrated 
in Fig. 5. 

The path integral of the axial force, along the atom tra- 
jectory in the optical funnel, reveals the potential barrier 
that needs to be overcome in order to reach the fiber. 
Consequently, all atoms funnelled toward the fiber expe- 
rience some deceleration that is dependent on their path 
trajectory. It is this complex interaction of their ini- 
tial state energies (transverse, axial and gravity poten- 
tial) with the varying forces of the optical dipolar funnel 
which determines whether they couple into the core. Fig. 
5 presents the axial velocities of each atom as a function 
of its distance from the fiber, in a Monte Carlo atom 
coupling simulation. The atoms trajectories have been 
labelled on this figure according to the three possible out- 
comes: escaped from the funnel (green), guided into the 
fiber (blue), or trapped in a dynamic balance between 
gravity and the axial dipolar force (red) . Once the guided 
atoms have gone through the repulsive axial potential 
barrier, they all start accelerating again inside the fiber 
by the action of gravity. The "trapped" atoms are pushed 
upwards and are recycled into the same near-triangular 
path of velocity-distance, being effectively trapped by the 
conservative potential. 



B. Atom dynamics in a microscopic optical funnel 

Fig. 6 displays five sequential snapshots of a 3D Monte 
Carlo simulation to fully illustrate the dynamics de- 
scribed. The figure is vertically compressed by a factor 
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FIG. 5. (Color online) The cixial velocity of each atom in 
the simulation is represented against its distance from the 
fiber. All atoms are released from 9 mm and exposed to a 
4 mW/300 GHz blue-detuned LGoi beam with a 20-/im waist 
radius within the fiber. Negative distances represent positions 
within the fiber core from the coupling input. Each velocity 
path is color-coded by the final outcome of its atom: guided 
(blue), trapped (red), and escaped (green). 



of 15 because the z axis spans over 17 mm while x and y 
axis are only 0.4 mm. In Fig. 6(a), the initial spatial dis- 
tribution of atoms is a 3D Gaussian function (spherical); 
although they initially appear to be in a disk distribu- 
tion due to the compressed z scale aforementioned. The 
confinement area is illustrated by the red cone colored 
from light yellow to red following the increasing trap po- 
tential. In our simulations, we observed that the mini- 
mal efficient coupling distance requires the radius of the 
maximum intensity of the beam, at MOT release loca- 
tion, to reach 1.8 times the 1/e-fold width of the MOT 
[as is the case in Fig. 6(a)]. This effectively corresponds 
to releasing from the MOT >98% of the atoms within 
the confinement limits of the funnel. 

Atoms initially outside the red cone will be immedi- 
ately lost. Those atoms inside the cone can escape only 
if they possess a transverse energy Exr greater than the 
trap potential depth, Uo{z). The transverse energy of an 
atom at position x, y, z is defined: 



ETr{x,y,z) = EKtivx.Vy) + U{x,y,z) 



(7) 



where Ext is the kinetic energy in the plane perpendic- 
ular to the fiber/beam axis and U{x, y, z) is the opti- 
cal trap potential at the atom position calculated from 
Eq (2). 

After the initial loss of the higher energy atoms in the 
distribution, in Fig. 6(b), the remaining cloud will no 
longer exhibit a Maxwell-Boltzmann distribution. We 
observe an areal compression of the cloud as it falls in the 
funnel. Although the total energy of each atom is con- 
served as it falls, the conversion of gravitational potential 
energy into kinetic energy leads to a large increase in the 
speed of the atoms. The diffraction of the funnel causes 




FIG. 6. (Color online) The dynamics of atomic coupling are represented with 5 snapshots taken during the simulation. The 
vertical axis is compressed by a factor of 15. The limits of confinement are represented by the outer light cone. The inside cone 
is delimited by the dark radius (peak intensity/e^). (a) Atoms at release time from the cooling lasers, 9 mm away from the 
fiber. Their density distribution is a 3D Gaussian but the compressed vertical scale makes it appear like a disk, (b) 16% are 
loaded in the hollow-beam funnel. Their phase space is compressed as they approach the fiber, (c) Almost all atoms approach 
the fiber within 0.5 mm, with some of them already coupling into the core, (d) Atoms are separated into two energy-selective 
populations: the coldest 8% are guided into the fiber core, the other 8% are pushed away by the axial component of the 
dipolar force, (e) Atoms relax up until gravity forces them down again toward the fiber. This cloud of atoms is trapped in the 
gravito-optical bottle trap as long as their total energy is conserved in the system. 



a conversion of the gain in axial veloeity into an increase 
in the transverse energy. The simulation shows that the 
increase in the transverse energy is proportional to the 
decrease in the area defined by the locus of points of max- 
imum intensity. That is, for a typical situation, with the 
MOT positioned 9 mm of the fiber, we find the radius of 
maximum confinement decreases by 160 times from the 
MOT release to the entrance of the fiber. The transverse 
energy of each atom also increases by 160 times. For- 
tunately, the trapping potential increases proportionally 
to the square of the decreasing fiber distance, which ex- 
actly matches the dependence of the transverse energy 
keeping the atoms in a trap of the same relative depth 
(ie. max transverse energy /optical potential depth is ap- 
proximately constant). 

Fig. 6(c) shows the atoms falling toward the fiber in a 
cone-shaped cloud, coming very close to the fiber surface, 
until a separation into two populations occurs [Fig. 6(d) 
and (e)]. The coldest fraction of atoms couples into the 
core; while others are repelled from the fiber. Continu- 
ation of the simulation shows that these atoms are exe- 
cuting the close-cycle trajectories shown on Fig. 5 as the 
"trapped" atoms. 

The principal determinant in the final state of the 
atoms is the initial transverse energy of the particle. The 
initial transverse energy of 1700 ®^Rb atoms was com- 
pared against their outcome in Fig. 7. It is shown that the 
coldest atoms (with lower initial transverse energy) cou- 
ple into the fiber core (blue circles). Atoms with higher 
energy immediately escape the hollow-beam at the MOT 
release distance (green triangles). Atoms with transverse 
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FIG. 7. (Color online) Initial transverse energy of 1700 Rb 
atoms labelled in function of their simulation results: guided 
into the fiber (blue dots), trapped by the gravito-optical bot- 
tle effect (red squares) or escaped the hollow-beam (green tri- 
angles). The 3D Monte Carlo simulations had the following 
parameters: 25 /jK, 200 fim diameter MOT launched at 9 mm 
into a LGoi beam of 40 fj.m diameter, 3mW and 300 GHz de- 
tuning. Only escaped atoms with transverse energy below 
8pK are shown. 



energy intermediate to those values penetrate the con- 
fining potential more deeply (red squares). These atoms 
have their axial velocity reversed before they reach the 
fiber. These atoms end up being trapped in an energy- 



selective gravito-optical bottle trap. 

Similar gravito-optical traps based on blue-detuned 
hollow-beams have been previously reported [39, 40], al- 
though their intention was to trap the atoms whereas we 
arc trying to avoid this. Under our set of conditions, 
the axial component of the blue-detuned repulsive force, 
opposing gravity, is only effective for atoms with higher 
transverse energy. During the coupling experiment, the 
gravito-optical bottle trap effectively behaves as a hot- 
atom filter. As the potential is increased in strength a 
larger fraction of atoms fall into this "trapped" class. Ob- 
viously, this bottle trap effect should be carefully man- 
aged in order to achieve optimal coupling of atoms into 
the fiber core. 

One can see from Fig. 7, that the transverse energy is 
not the singular determinant of the eventual classifica- 
tion of the atom. Wc find that the initial axial velocity 
can have a minor influence for those atoms with axial 
velocities at the cxtrcma of the distribution. This causes 
an atom to fall into a class which was unexpected from 
their initial transverse energy alone. In Sec. IV, we will 
consider how the experimental parameters influence the 
atom's distribution among the three classes. 
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FIG. 8. (Color online) Atom coupling efficiency as a function 
of tiie force parameter of the blue-detuned hollow-beam and 
the distance the atoms in the MOT are released from the fiber. 
The optimal line (blue circles) shows the force relation with 
MOT distance to realize the best coupling efficiency. The low- 
force line (green squares) and high-force line (red triangles) 
represent the force-distance combination that would deliver a 
factor 2 reduction in the coupling efficiency. 



IV. SYSTEM PROPERTIES: OPTIMIZING 
COUPLING EFFICIENCIES 

In this section, wc will focus on evaluating the coupling 
efficiency of our specific arrangement. In particular, we 
will discuss the dependence of the efficiency on: the light 
force parameter, MOT-fiber distance and fiber core di- 
ameter. We will also make a comparison of this hollow- 
core approach with an optimal coupling scheme based on 
rcd-detuncd Gaussian beam coupling. 



A. Optimal coupling efficiency discussion 

In our simulations with ^^Rb, we selected realistic 
MOT parameters of T=25 ^K and the trap radius a = 
70/im using a standard setup [35]. We have centred the 
MOT on the fiber axis and made use of a HC-PCF that 
guides a 40-/im diameter LGqi beam. The remaining 
degrees of freedom are the dipolar light force parameter 
(k = Pj IS) and the distance between the fiber and the 
centre of the MOT. The simulation shows that there is 
an optimal potential for maximum coupling at a given 
fiber-MOT spacing. This optimal value is shown by the 
blue circles line on Fig. 8. Wc also show two other lines 
in which the coupling fraction has fallen by a factor of 
two from the optimal value. The relationship between 
the force necessary to maintain optimal coupling and the 
MOT distance is linear. This arises because we need to 
compensate for the additional diffraction by increasing 
the beam power to maintain the equivalent initial con- 
ditions. The strength of the axial potential barrier that 
increases proportionally with k, is a determinant factor 
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FIG. 9. (Color Online) Monte Carlo simulation of the optimal 
coupling efliciency of a 25 /iK, 200 /im diameter **^Rb MOT 
launched from a fixed distance into: a 40-/im diameter (blue 
circles) and 50-/im diameter HC-PCF (red squares) guiding a 
blue-detuned LGoi beam. 



in the coupling efficiency. Its effect on the coupling effi- 
ciency is compensated by the linear increase of gravita- 
tional potential energy with distance. The axial barrier 
stays therefore at the same relative strength for a given 
atom and the optimal coupling efficiency is maintained. 



1. Optimal coupling as a function of the MOT-fiber 
distance 

Fig. 9 shows the maximum atom coupling efficiency 
obtainable as a function of the distance from the MOT 
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FIG. 10. (Color Online) Releasing the MOT at 9 mm from 
the fiber, the coupling efficiency of *^Rb atoms is illustrated 
in relation to the light force parameter (blue circles). The 
trapping efficiency is also represented against the same light 
force parameter (red triangles). Best coupling efficiency is 
obtained with the maximum force allowed before increasing 
the fraction of atoms trapped in the gravito-optical bottle 
effect 
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FIG. 11. (Color Online) On the left axis, minimum coupling 
distance to achieve the best coupling efficiency is represented 
as a function of the fiber core radius (red triangles) . On the 
right axis, the best coupling efficiency achieved with a partic- 
ular fiber core radius is shown (green circles) . The red dashed 
line is a linear fit of the distance-radius relationship and the 
green dashed line is a quadratic fit to the coupling efficiency. 



for two different fiber core radii. The analysis reveals, 
that beyond a minimum distance, the best coupling ef- 
ficiency can be achieved. This value can be maintained 
for increasing distance. This distance is set when the 
radius of the diffracting funnel is at least 1.8 times the 
size of the MOT. Below this threshold the Gaussian dis- 
tribution of atoms in the MOT is trimmed by the funnel 
confinement limits. Even if the coupling efficiency can be 
maintained with longer distance, it is not desirable to do 
so because the increased light force parameter will results 
in higher scattering rates in the fiber core. Therefore, it 
is best to work with the lowest force possible by selecting 
the shortest distance for optimal coupling efficiency. The 
minimum distance consistent with capturing the maxi- 
mum fraction of the MOT is 9 mm with a 40- /im core 
diameter and 11.25 mm with 50- /xm. 



2. Coupling efficiency as a function of the light force 
parameter 

The fraction of in-coupled atoms as a function of the 
light force parameter in the hollow-core guide is displayed 
on Fig. 10 for our standard set of conditions. We have 
chosen a fixed distance of 9 mm between the MOT and 
the 40-/im fiber. The coupling efficiency increases up 
to a certain force, and then, decreases as the fraction of 
trapped atoms grows. We see that there is an optimal 
depth to the potential: this characteristic appears be- 
cause we need to set a balance between the depth of the 
potential, in terms of atom capture, against the appear- 
ance of a population of trapped atoms. The coupling- 
efficiency-light-force relationship shows the same opti- 



mum with different MOT release distances and different 
fiber core radii. Thus, for a set distance and fiber core 
radius, the optimal coupling efficiency is always achieved 
with the maximal force just before the appearance of 
trapped atoms. The fraction of trapped atoms rises much 
faster because not only guided atoms are converted into 
trapped atoms but the increase in force also increases the 
fraction of captured atoms. Nonetheless, their higher 
transversal energy means they cannot be converted to 
coupled atoms either. 



3. Scaling of coupling efficiency with fiber core diameter 

If the objective is to couple large numbers of atoms into 
a fiber, it is clearly an advantage to work with larger core 
fibers. With a larger core, there is less diffraction (diffrac- 
tion angle is inversely proportional to the radius) and this 
results in a naturally weaker axial force. Comparing dif- 
ferent fiber core radii, the maximum coupling efficiency 
and minimum coupling distance to achieve that efficient 
coupling are displayed for a range of realistic fiber core 
sizes in Fig. 11. Our simulations have shown that the 
maximum coupling efficiency is linearly proportional to 
the fiber hollow-core area (a quadratic relation with core 
radius). The minimum coupling distance is also a lin- 
ear relation of the core size. For best coupling efficiency, 
we previously observed that the minimum coupling dis- 
tance is achieved (with a LGqi beam) when the radius 
of maximum intensity of the diffracting hollow-beam is 
equal to 1.8 times the diameter of the MOT. Hence, for 
this fixed funnel radius, the relationship of the distance 
with the fiber core radius is d = 1.8cr/tan{d), where d 



is the coupling distance and 9 is the diffraction angle of 
the hollow-beam (inversely proportional to the core ra- 
dius). The minimum distance is therefore a linear func- 
tion of the fiber core. With increasing core dimensions, 
the maximum light force that can be applied is limited by 
the axial dipolar force (the gravito-optical bottle effect), 
as reported in Sec. IV A 2. It is observed that increasing 
K proportionally to 7^^ will maximize the coupling effi- 
ciency. The quadratic increase in the coupling efficiency 
with core radius (shown in Fig. 11) is a combination of 
the decreasing axial force by 1/r and the linear increase 
of the gravity potential energy (due to increasing linearly 
the minimal distance with r). 



B. Optical scattering in the guiding beam 

We consider the difference in light scattering between 
atoms guided in a conventional rcd-dctuned Gaussian 
beam and a bluc-dctuncd hollow-beam. For frequency 
a detuning large compared to the atomic linewidth and 
at low saturations, the scattering rate per atom can be 
expressed: 






(8) 



In both guiding arrangements, the cold atoms will 
be concentrated around the central axis depending on 
their initial transverse energy. In the blue-detuned guide, 
the atoms will spend the predominant time in the dark 
whereas, in the Gaussian beam, they are exposed to the 
maximum intensity. We have compared the light scatter- 
ing rates of atoms in repulsive and attractive potentials 
during the process of coupling and then, guiding inside a 
HC-PCF. We have used an identical maximum potential 
and a detuning of 300 GHz from resonance in both cases. 
We have calculated these light scattering rates using two 
different methods. 

In the ffist method we explicitly calculate the local in- 
tensity, and hence scattering rates, at each time step of 
our modelled trajectories in the Monte-Carlo simulation, 
while in the second we have calculated the average light 
intensity in the volumes explored by atoms of a specific 
energy and then used this to represent the average scat- 
tering rates. Both techniques predict that only a small 
fraction of scattering events (less than 8% of the total) 
occur during the coupling phase and hence we will ig- 
nore this phase in our calculations. Fig. 12 compares the 
relationship between light scattering rates and the cou- 
pled atom transverse energy, using this first technique, 
for a Gaussian (black squares) and a LGqI trap beam 
(green triangles). We observe the expected behavior that 
a hollow-beam guide has much lower scattering rates than 
a Gaussian guide except for atoms with a transverse en- 
ergy close to the guiding potential. 

On Fig. 12 we also display the expected scattering rates 
from a simple model that averages the beam intensity 
over the area of confinement of atoms of a particular 
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FIG. 12. (Color Online) Light scattering rates of cold atoms 
guided in a HC-PCF by an attractive Gaussian potential (up- 
per red curve and black points) and by a repulsive hollow-core 
potential (lower blue curve and green points). We have cal- 
culated the expected scattering rates by a statistical method 
(smooth curves) and by examination of the local intensity 
through the explicit trajectory of the atoms (points). The 
two methods show excellent agreement. Simulated param- 
eters are: 9-mm MOT-fiber distance, 40-/im core diameter, 
0.57 mK trapping potential depth in the core, and 300 GHz 
detuning. 



transverse energy. The radius of that area is determined 
by a trapping potential that is equal to the transverse en- 
ergy of the atoms. The calculation of the light scattering 
rates for the Gaussian (red line) and LGqI beam (blue 
line) shows excellent agreement with the light scattering 
rates by the explicit technique. Thus, this investigation 
shows that blue-detuned hollow-beam guides arc supe- 
rior in the guidance of cold atoms because of the reduced 
scattering rates while maintaining equivalent guiding po- 
tential. It was also observed that for different optical 
trap depths, the ratio of scattering rates between Gaus- 
sian and hollow-beam, for a given transversal energy, is 
maintained. 



Model validation and comparison with 
red-detuned Gaussian beams 



We tested our model against the experimental data of 
atom guiding efficiency from Mestre et al. [27]. The 
Rb atoms were guided using collimated blue-detuned LG 
modes of order 1 to 12. Because their model didn't take 
into account light scattering losses, they added an em- 
pirical loss factor fitted against their experimental data. 
We applied the same light scattering loss factor to our 
simulated results. Table I shows that we obtained ex- 
cellent agreement with the experimental data. We also 
applied our model against the cold atom coupling experi- 
ment reported recently by Bajcsy et al. [25] (summarized 
in Sec. I). We modelled their initial loading configuration 
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Beam order Exp. guiding eff.(%) Our model calculation (%) 

O 
8.1 
12.75 
18.3 
19.1 



TABLE I. Comparison between experimental data from 
Mestre et al. [27] and our model prediction in the guiding 
efficiency of rubidium atoms in coUimated blue-detuned LG 
beams of various orders. 



as a spatially uniform disk distribution of 440 yuni diam- 
eter, 1 mm from the fiber, with a temperature distribu- 
tion equal to that of the MOT. We obtained a coupling 
efficiency of 1.5%, with our model compared with their 
0.3% experimentally. Our calculation has not assumed 
any losses in the MOT loading process and no light scat- 
tering losses during the red-detuned Gaussian coupling. 

When comparing the coupling efficiency of the blue- 
detuned hollow-beam with a red-detuned Gaussian 
beam, the latter has the advantage of a larger confin- 
ing area (no confinement limits) and the axial dipolar 
force is attractive, which facilitates coupling. Using the 
same fiber and identical trap depth potential in the guide, 
the best coupling efficiency obtained with a red-detuned 
Gaussian beam has shown to be a factor 4 more efficient 
(46%) than the blue-detuned LGqi beam (11%) when us- 
ing a 40-/im diameter fiber. 

Nonetheless, the blue-detuned hollow-beam is poten- 
tially more useful in some circumstances for guidance of 
cold atoms because of reduced light scattering rates for 
an identical trapping depth (see Sec. IV B). 



V. CONCLUSION 

In this paper, we presented a model and simulations 
aimed at efficient coupling of atoms from a MOT into 
a hollow-core photonic-crystal fiber using a blue-detuned 
LGoi beam. In contrast with previous approaches, only a 
single beam is required to achieve both efficient coupling 
and guiding within the fiber. The key goal is that the 



atoms are guided in locations with minimal light intensity 
in order to minimize shifts of the internal states and from 
decoherence associated with light scattering. 

We showed that, with a blue-detuned hollow-beam fun- 
nel, there is a minimum efficient coupling distance, which 
can be simply calculated from the measured parameters 
of the MOT and the fiber core radius. Our model iden- 
tified a new gravito-optical bottle trapping effect that is 
significant for the coupling dynamics, when using a low- 
order, blue-detuned hollow-beam funnel. This trapping 
effect must be carefully managed to optimize coupling 
efficiency. We identified that optimal atom coupling con- 
ditions, for a fixed MOT-fiber separation and fiber core 
radius, is achieved when using the highest light force pos- 
sible before the appearance of trapped atoms. This result 
disputes the naive expectation that strong guiding forces 
will lead to maximum in-coupling to the fiber. In con- 
tradiction with this expectation, we find that there is an 
optimal potential depth that will maximize the coupled 
fraction from the MOT. 

When taking into account experimental constraints: a 
minimum fiber-MOT separation of 9 mm, hollow-beam 
power <50mW; we obtained optimal coupling efficiency 
of 11% with a LGoi beam when using a 40-^m core diam- 
eter, 25 mW of power and far-detuning of 3THz. Larger 
core fibers could achieve better coupling efficiency but 
would require higher light force and experimental con- 
straints would force trading-off power for a reduced de- 
tuning, which is not desirable to minimize light scattering 
losses. 
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